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PREFACE 


The work described in this report was performed by the Applied 
Mechanics Division of the Jet Propulsion Laboratory, as part of the NASA 
support of the Helios Project* 
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ABSTRACT 


This report presents methods of measuring moments of inertia with 
very high accuracy. 

The moment of inertia of the Helios Spacecraft about its spin axis was 
determined by use of a Tr roll-fixture n using two sets of crossed flexure 
pivots as elastic constraints. The test procedure entailed measurement of a 
system oscillation period with each of a set of added moment-of-inertia 
increments. The tare effect of the fixture was determined by a like process 
and was subtracted from the gross value to yield the spacecraft roll 
moment of inertia to an estimated accuracy of 0. 2%. 

"Lateral" moments of inertia (i. e. , about each of three axes normal to 
the spin axis) were determined by a gravity pendulum method that makes 
use of the fact that any physical pendulum has a minimum period of oscilla- 
tion determined by a particular distance from the axis of rotation to the 
system center of gravity. In situations where a knife-edge support is used, 
this distance is equal to the system centroidal radius of gyration. In the 
subject tests, the pivoting action was provided by hardened pins rolling on 
flat ways. The effect of the finite radius of the pins was considered in 
deriving the equations of motion, from which an error analysis revealed the 
criterion for maximum accuracy in determining the square of the centroidal 
radius of gyration. 

The swing fixture provided for a number of optional pivot-pin locations 
giving precisely known distances between successive axes of oscillation. 

This fixture, with provisions to support the spacecraft, was ballasted to 
bring its vertical c.g. close to that predicted for the spacecraft. This 
ballasting was done not because the test method requires an accurate 
foreknowledge of specimen c.g. position but, rather, to minimize errors 
in the par allel- axis transfer term while removing the tare of the fixture. 

Though the centroidal moment of inertia of the swing fixture was over 
twice that of the spacecraft, an error analysis showed that accuracies of 
better than 1.0% were realized for the two lateral principal moments of 
inertia. 
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SUMMARY 


Principal moments of inertia of the Helios Prototype Spacecraft have 
been determined at JPL from measurements made during the period, 

July 2 through July 8, 1974. 


Measurements were made with the magnetometer booms installed and 
in their launch configuration. For measurements of moment of inertia about 
the roll (spin) axis, the antenna reflector was taped to the spacecraft to lock 
out M despin M freedom. 

Following are the mass properties determined in this program, with 
estimates of probable errors: 


W = 349. 3 kg (±0. 1% max. ) 

Z = 929 ±19 mm (above separation plane) 
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This report describes the program of test and analysis leading to these 
results . 


INTRODUCTION 

This report documents the determination of the principal moments of 
inertia of the Helios Prototype Spacecraft. 

At the outset of the program the following requirements were defined: 

An accuracy of \% in spacecraft moments of inertia was 
stipulated . 
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A maximum period of six weeks was available for test planning, 
fixture design and fabrication, and testing. 

The determinations were to be independent of prior estimates of 
spacecraft inertial properties. 

The spacecraft was required to remain vertical, or nearly so. 

A consideration of all of these factors, coupled with a review of 
spacecraft nominal inertial properties as presented in Table 1, led to the 
choice of a pendulum method for determining radii of gyration about axes in 
the X-Y centroidal plane and a method using elastic restraint for the deter- 
mination of the moment of inertia about the Z (roll) axis. 

The particular pendulum method selected had been suggested in 
earlier analytical work as one not requiring prior knowledge or independent 
measurement of the vertical location of the center of gravity and one having 
a high accuracy potential. However, it had not been assessed experimentally 
and so, in conjunction with the primary effort, a scale model program was 
instituted, as reported in Appendix C. 

Common to both methods employed in these tests is the requirement 
for accurate measurements of mass (M), length (L) and time (T). For the 
pendulum- method tests, weights were determined by use of precision- 
calibrated load cells with read-out accuracies to 0,05 kg and probable 
errors not exceeding 0. 02% of full scale. The weight increments used in 
the roll moment-of -inertia tests were measured on a certified gram 
balance to the nearest 0. 1 gram. 

Critical lengths were measured to an accuracy of 0. 06%. 

Elapsed times for a predetermined number of cycles of oscillation 
were measured with a photoelectric timing system giving six significant 
figure read-out. 

The influence of these accuracies on the overall accuracy determina- 
tions is reported in Appendices A and B. 

The following sections of this report deal with the theory, the 
implementation, the results and the assessments of the determinations of 
the principal moments of inertia of the Helios Prototype Spacecraft. 
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DETERMINATION OF THE INERTIA ELLIPSE IN THE X-Y CENTROIDAL 
PLANE 


The Theory 
a) Equation of motion 

A pendulum method was used to determine moments of inertia of the 
spacecraft about axes in the X-Y centroidal plane, 

Kater’s reversible pendulum has been classically used in the deter- 
mination of the gravitational constant to six significant figures, showing the 
extreme accuracy obtainable in a pendulum experiment. 

The ideal physical pendulum is characterized by its frequency 
equation. 


2 

u) 



( 1 . 1 ) 


where 

circular frequency, radians/sec. 
gravitational constant, m/sec , 

distance, in meters, from the center of rotation to the "center 
of oscillation, 11 also referred to as the n center of percussion. IT 

The parameter, D, is recognized as the length, L, of an equivalent 
simple pendulum. It is given by 


w = 

g = 
D = 


D 


2 2 
r + p 

o 


(1.2) 


where is the centroidal radius of gyration and r is the distance from 

the pivot point to the center of gravity (c.g. )• D has a minimum and, 

consequently, the frequency of oscillation has a maximum for r = p^. 

An error analysis shows that the maximum accuracy in the determination of 

p ^ exists also for r = p . 
r o r o 
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In a test situation axes of oscillation are commonly realized with 
knife-edge pivots. In dealing with massive systems, the knife-edge design 
poses the materials problems of galling and fretting under bearing stresses. 
Accordingly, an alternative is used in which the pivot system is comprised of 
finite-diameter pins rolling on flat ways. While this system provides a 
reduction in bearing stresses, it also poses some readjustments in the 
equations of motion, as the n rest point' 1 now describes cycloidal motions 
under angular oscillations. 

Figure 2 shows the geometry of such a system, in which 
a = radius of pin 

r = distance from pin support to the system c. g. 

6 = angular amplitude of displacement 

x = horizontal displacement of system c.g. 
z = vertical displacement of system c.g. 

m = system mass 

= system centroidal radius of gyration 

For a small angular displacement, 0 (assuming no slippage of the 
pin on its support), the horizontal displacement of the c.g. is 

x = (r + a)0 - a0 - a( 0 - sin 0 ) (1.4) 

or, to a first approximation, 

e 3 

x « re - a-|- (1.5) 


The associated vertical displacement is 

z = (r + a)( 1 - cos 0 ) 

2 

~ (r + a)-^- 


(1.6) 
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The kinetic energy of this system is 


„ 1 21 2-2 
T = mv + m P Q 6 


(1.7) 


where 


2 * 2*2 
V = X + z 


( 1 . 8 ) 


Use of Eqs. 1.5, 1.6, and 1. 8 in Eq. 1 . 7 leads to 


1 tA2 


T = -g- 16 4- higher order terms in 


(1.9) 


where 


I = m (r 2 + p Q 2 ) 


(1. 10) 


Neglecting the higher order terms leads to the inertia moment, 


d /9T\ 
dt \d¥] 


ie 


(i. ii) 


The potential energy of the system is (for small amplitudes) 


V = mgz 


0 

~ mg {r + a) 


( 1 . 12 ) 


The restoring moment is 


9V 

80 


mg (r + a)6 


( 1 . 13 ) 
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The equation of motion without external forcing and neglecting damping 


is : 


16 + mg (r + a)0 = 0 


For simple harmonic motion, 


Then 


or 


0 = 0 e 

o 


iwt 


Thus 


Let 


0 = -oo 0 



D = 


2 , 2 
r + p 
o 

r + a 


w 


£ 

D 


T 2 = D 


(1.14) 


(1.15) 


(1.16) 


d.l) 


where T is the period of oscillation in seconds* 
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It is significant to note that the pin radius, a, enters only into the 
potential energy relation of Eq. (1. 12); it does not appear in the kinetic 
energy relations (Eqs. 1.9 and 1. 10) because the instantaneous center of 
rotation remains at the interface between pin and support. 

2 

The method used to determine p^ was to perform measurements of the 
period of oscillation for several locations of the axis of oscillation (i. e. , 
several values of r) in a search for the minimum period, T^, and, corre- 
spondingly, a minimum value of D. 

2 

b) Determination of p Q from a minimum of D 
In Eq. (1, 16), let r + a = R, so that 

r = R- a (1.17) 

Substituting of Eq. (1. 17) into Eq. (1. 16) leads to 


D 


R - 2a + 



( 1 . 18 ) 


For the minimum sensitivity of D with respect to R, i. e. , 


dD 

dR 


0 



Substitution of Eq. (1.19) into Eq. (1. 18) leads to a minimum value, 
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from which 



( 1 . 20 ) 


c) Determination of from a best-fit curve of period vs pivot pin 
location . 

2 

The square of the radius of gyration, p Q , can also be determined 
from a best fit, in a least-squares sense, of the test data, i. e. , the oscil- 
lation periods, T., as a function of the oscillation axis locations, r.. 


P 


o 


2 


n 


n 



+ a ) — ^2 
4tt 




( 1 . 21 ) 


where n is the number of axis locations as shown in Appendix A. 
2 

d ) Error in from a minimum in period of oscillation 

Differentiating Eq. (1.20) gives 


d p 
r o 


2 


P 

o 


2 




da 

a 


Eq. (1.1) may be written in the form. 


D 


- — T 2 


m . 2 m 
4tt 


( 1 . 22 ) 


where is the minimum period of oscillation. 


dD m 1 / aD m\ , , 1 / 9D m\ _ 

D m \ / dg D„ (“8X - ) dT m 
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Thus 
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(1.23) 


In the special case of the knife-edge pivot (a = 0), Eq. (1.23) reduces 
to 


A 



P 


2 

o 



AT 


4- 4 


m 


m 


(1.24) 


The error in p is not explicitly related to the radial distance, r. 
° 2 

This makes the determination of independent of the center of gravity 
to the first order of approximation. 


Once has been obtained, the radial distance, r, from the system 
c.g. to the pivot point associated with may be determined from 
Eqs. (1. 17) and (1. 19)- 
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e) Error on p Q from a best-fit curve of period vs. pivot pin location 

Appendix A presents a statistical error analysis based on the best-fit 
curve. An estimate for the mean square error is (Appendix A, Eq. A- 18) 


S 


2 





- <r, + a) 




2 

+ r . 

i 


2 


Test Implementation 
a) Setup 

A support cradle accommodating a spacecraft adapter is shown in 
Fig. 3. The box section supporting the adapter was designed for high tor- 
sional stiffness; even though the inertia torques imposed during oscillation 
were calculated to be quite low, the associated angular deflections between 
sidearms were held to very small values (less than 0. 00001 rad. ) to limit 
periodic vertical migration of the instantaneous center of rotation to less 
than 1 mm with peak amplitudes of oscillation less than 1 deg* 

Likewise, the sidearms carrying the support pins were designed for 
high stiffness in the plane of oscillation. 

The bare cradle was provided with ballast weights, based on mass 
property calculations, to raise the center of gravity of the cradle assembly 
(plus adapter and V-band) close to the nominal vertical c.g. of the space- 
craft in order to minimize the transfer term as seen from Eq. (1. 27). 
These ballast weights also provided a series of precisely spaced holes into 
which the pivot pins could be inserted. (See Fig. 4..) 

Hardened and ground bearing plates, bolted to support stands, pro- 
vided for rolling contact in the cradle suspension. 

A specially-designed hoisting sling suspended from a bridge crane 
was used to lift the entire assembly during the operations of weighing and 
of changing pin positions. 

Prior to installing the adapter on the cradle, the test system was 
proof-loaded with a dead weight of 1360 kg. 


10 
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After installation of the adapter and V-band, the carriage assembly 
was weighed by use of a precisely calibrated load cell between the crane 
hook and the hoisting sling. The tare weight of the hoisting sling was also 
recorded. 

The vertical position of the c, g. of the cradle assembly was deter- 
mined by measurement of the static moment required to rotate the assembly 
through 90 deg. Accuracy of the angle of rotation was determined by a 
Hilger and Watts inclinometer giving a vernier reading to the nearest min- 
ute of arc. 

The data obtained from these static measurements of cradle weight 
and c,g. determination are reported below, under Test Results. An error 
analysis is presented in Appendix A, 

Mechanical stops were provided to prevent cradle rotation during 
installation of the spacecraft. Following spacecraft installation, the weight 
of the entire assembly was determined. Views of the Helios spacecraft on 
the test fixture are shown in Figs. 5 and 6. 

b) Tests 

There were four series of oscillation tests conducted. One of these 
was to determine the centroidal radius of gyration of the fixture alone. The 
remaining three included the spacecraft with oscillations about the x, y, 
and p axes, respectively. Each series was repeated at least once to 
evaluate reproducibility. 

The following features were common to all of these tests. The photo- 
electric sensor was positioned under one end of the cradle platform, which 
carried a plate that interrupted the light beam twice per oscillation (Fig, 7), 
With the cradle at rest, the sensor was placed so that the electronic counter/ 
timer (Fig. 8) was on the verge of tripping, (The resolution in manual place- 
ment is estimated to be about 1,5 mm. The associated error ascribable to 
damping of amplitude over the elapsed time for 40 cycles of oscillation is 
estimated to be less than 0. 002%. ) 

In each test series, timings were made with five different (parallel) 
axes of rotation by use of five consecutive hole positions for the pivot pins, 

*The n p n axis is rotated 33. 75° from the -fy axis toward the +x axis. 
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The positioning of these pin locations in relation to the spacecraft separation 
plane is shown in Fig, 9. With each pivot axis, two measurements were 
made of the elapsed time (for 40 cycles of oscillation) with an initial ampli- 
tude at the sensor of about 25 mm, and another two were made with an 
initial amplitude of about 13 mm. An amplitude effect on elapsed time was 
discernable and the times for Tl zero amplitude" were obtained by graphical 
extrapolation. In this way external aerodynamic effects and small non- 
linearities in kinematic behavior have been essentially eliminated. 

One deficiency of the fixturing is to be noted. Each bearing plate 
on which the support pins rolled was canted because of warpage in welding 
its faceplate to the central column. Pressure of time precluded correcting 
this situation by a final machining, and the result was that each pin did not 
have line contact, but, rather, bore on one edge. This resulted in pro- 
gressive galling of the pins, which subsequent Rockwell tests have shown 
were not as hard as the bearing plates. In future tests of a similar nature, 
this problem would be eliminated. 

Test Results 

a) Configuration 1, cradle plus adapter and V-band ; 

Weight (including hoisting sling) 469. 7 kg 

Tare of hoisting sling 51,9 kg 

Net, cradle plus adapter and V-band 417. 7 kg 

2 

C.G, determination 



The V-band was not available at the time of these measurements. Its 
effects are accounted for in the following calculations. 
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= 186.7 kg (less V-band) 

W i = 417.7 kg (including V-band) 

L = 2. 169 m 

AW = 7. 54 kg @L r = 2. 040 m (V-band) 

The static moment about the center of pin D is 

M = W. , r ' - (PL + AWL’) = 0 
tot o 

PL + AWL 1 


r 1 = 1, 0064 m 

o 

From the bas e of pin B, 

r = r ' - 0. 1022 
o o 


r = 0. 9042 m 
o 



Weight (including hoisting sling) 820. 1 kg 

Tare of hoisting sling + load cell cable 53. 0 

Net, cradle plus spacecraft 767. 1 kg 

Cradle plus adapter and V-band 417.7 

Net, Spacecraft 349. 3 kg 


c) Timing measurements, all configurations 

The results of the timing measurements are presented in Tables 2 
through 5. 


Data Processing 

By gravimetric measurements conducted by the Caltech Seismological 
Laboratory, the gravitational constant on the ground floor of JPL's Gyro 
Laboratory has been established as 

2 

g = 9.79543 m/sec 
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In consideration of the small elevation difference between the Gyro Lab 
and the site of the subject tests (Bldg. 144), the value listed above has been 
used, rounded off to five significant figures. 

Appendix A makes use of the test data presented in Tables 2 through 5 
to obtain, through an iterative process, the "least squares" best fit of dis- 
crete data points to the theoretical relationships. Here, r is replaced by 
r Q + , where, as used, r Q is the distance from the base of the pivot pin 

at Hole B to the system vertical c.g. position. This process requires 

2 

initial estimates of r Q and Pq for each data set. These initial estimates 

were obtained from plots of oscillation period, T vs x. The minimum time, 

T , obtained for each configuration from a visually-faired curve was used 
« 2 

in Eqs. (1. 22) and (1. 20) to obtain a trial value of p . The value of x at 

o 

which this minimum appeared was used to obtain the trial value of r . 

o 

The initial statistical curve fitting reported in Appendix A has been 
used primarily to obtain improved values of r . These, of course, should 
be identical, regardless of the spacecraft orientation on the cradle. How- 
ever, there were differences in the three values of r Q with the spacecraft 
x, y, and p axes parallel to the oscillation axis, although the maximum dif- 
ference from the average was only 0. 7%. A "consistent" average value of 
r for the three axes was chosen for the final calculations of p for each 

U r o 

test value of x.. These values were used in a least- square s best fit, 

1 2 

Eq. (1. 21), to the linear relationship, p = a constant, for each spacecraft 

orientation. Additionally, the variance has been calculated to obtain a most 

probable error in each value of p 

r o ' 

The results of the processing of the data for the fixture only and for 

the fixture plus spacecraft in each of the three orientations are presented 

in Table 6. Results of the final calculations are presented graphically in 

Figs. 10a through 1 Od as plots of T vs. x. Here, the test points are circled 

and the theoretical best fits are shown as solid lines. The results are also 

shown in Figs. 11a through lid wherein the dashed lines represent the 

2 

statistically derived values of p , showing that the dispersion from the 
mean is random with, as expected, no apparent correlation between p 

° 

and the location of the axis of oscillation. 
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a) Separation of spacecraft and fixture 


To obtain the desired mass properties of the spacecraft (subscript S), 
the effects of the fixture (subscript F) must be removed from the total 
(subscript T). 

Letting M, r^, and represent mass, c. g. coordinate, and radius 
of gyration, respectively, there exists, for static equilibrium, the relation 


from which 


M F ^oT " r oF^ + M s ' r oT " r oS) ° 


M, 


r oT " r oS " M„ Vi oT ” ‘oF' 


(r„rp - r . 


(1. 25) 


also, the total moment of inertia is 


(Mp + M g ) p oT Mp p Q p + Mp (r oT - r Q p) 


+ M S p oS + M S (r oT " r oS ) 


(1. 26) 


From Eqs. ( 1 . 25) and ( 1 . 26 ) , there follows , 


2 

P oS 


/ M \ M M 

\ l + m|) P oT-Mf P oF - < r oT - r oF> Vq 



(1. 27) 


Because of the fixture ballast weights, r ^ - r is very small (cf. 

o x or 

Table 6), and the last term may be neglected. 

2 

Calculations of the variances on are presented in Appendix A, 
together with moments of inertia about the x, y, and p axes and their 
associated probable errors. These results are listed in the Summary.. 
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b) The inertia ellipse 


R 



2 2 

I = I cos a + I sin a - TJ sin 2 a 
p x y xy 


(1. 28) 


where = product of inertia about the x and y axes. From Eq. (1.28), 


U 


xy 


2 2 

I cos a -f- I sin a - I 

_JE Y E 

sin 2 a 


(1. 29) 


Also, for other orthogonal axes, r and s, 


U = x (I - I ) sin 2 <j> -f- U cos 2<b 
rs 2 x y xy T 


(1. 30) 


where 4> is the angle between the r and x axes. For the r and s axes to be 
principal axes, = 0, and 


tan 2 4> 


-2U 



I - I 
x y 


(1.31) 


of 


Eq. 


The procedure, here, is straightforward. U is calculated by use 
(1.29) for obtaining <j> by Eq. (1.31), Then Equation (1.28) is used. 
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with 4> replacing or, to obtain the principal moments of inertia* An error 
analysis is given in Appendix A* 

Here, 

I = 181. 69 kg - m^ 

I = 186. 57 kg • m 2 

y 

I = 188. 27 kg • m 2 
P 

a = 56.25° (counterclockwise from +X axis) 

from which 4> can be calculated. 

<j> = 27. 445° = 0.4790 rad. 

Then from Eq. (1. 28), with <j> replacing a, 

I = 188. 37 kg . m 2 = I 

r & max 

and 

I + I - I = I . 
x y r min 

179. 88 kg • m 2 



These results are listed in the Summary. 


DETERMINATION OF THE MOMENT OF INERTIA ABOUT THE Z (ROLL) 
AXIS 


Measurement Concept 

For a body in free sinusoidal oscillation about a fixed axis, the period 
of oscillation, T q is determinable from the relation, 


w 


2 

o 


4tt 2 K 



o 


( 2 . 1 ) 
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where K is the restoring moment per radian angular displacement and I 
is the moment of inertia of the body. 

If an increment, AI of the moment of inertia 1^, is added to the system, 
the new period, T, is obtained from 


4tt 2 K 

T 2 = + AI 


( 2 . 2 ) 


Dividing Eq. (2. 1) by Eq. (2. 2) leads to 



(2.3) 


Thus the square of the period is a linear function of AI in which the 
slope mis: 


m 



Thus 


I 

o 



o 


m 


(2.4) 


(2.5) 


This approach, which has been in occasional use, entails merely the 
precise determination of several AI T s and the associated periods of oscil- 
lation from which to establish the slope, m, and thus 1^ from Eq. (2. 5). 

Test Implementation 

The easiest implementation of the concept is to use elastic restoring 
torque. Accordingly, a test fixture was designed to employ two vertically 
separated sets of crossed flexure pivots between a floor -mounted support 
stand and the rotational table. This fixture, shown in Fig. 12, carries a 
long transverse arm with provisions for attaching weight increments. 
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The flexures were designed to withstand the dead load of the spacecraft 
and the bending stresses due to maximum oscillation amplitudes with an 
ample margin of safety. The lower pair of flexures is shown in Fig. 13. 

Prior to installation of the spacecraft adapter on the rotational table, 
the system was proof-loaded with a dead weight of 136 0 kg, and oscillated 
through the design amplitude while under load. 

A test was also made to ascertain that the axis of table rotation was 
actually on the spacecraft roll axis as established by the adapter position. 

With the spacecraft adapter and V-band in place, tests were made to 
determine the moment of inertia of the fixture. . The spacecraft was then 
installed, and timing measurements of the complete system were made. 

Figs. 14 and 15 show the spacecraft in position on the test fixture. 

Fig. 16 shows the placement of the photoelectric timing system and 
four of the weight increments secured near the end of the transverse arm at 
a precisely determined distance from the roll axis. 

Test Results 

The data obtained in these tests are presented numerically in Tables 7 
and 8 and graphically in Figs. 17 and 18. 

Data Processing 

The test data of Tables 7 and 8 have been used in a statistical analysis 
presented in Appendix B to derive the roll moments of inertia of the fixture 
and the spacecraft plus fixture, and the statistical variances on each. The 
results are presented in Table 9* 

CONCLUSIONS 

The two test methods used for obtaining principal moments of inertia 
are deemed wholly satisfactory. The probable errors in spacecraft princi- 
pal moments of inertia are below the target values of one percent set at the 
beginning of the program in spite of the large tare moment of inertia that had 
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to be accepted in the pendulous tests. The test hardware provided sufficiently 
low damping (damping coefficient, y < 0. 01%) that errors in elapsed times 
of oscillation contributed negligibly to overall error. 

Although the results are excellent, the full accuracy potential of the 
multiaxis, pendulous method employed here was not completely realized 
because of alignment and material deficiencies local to the support pins and 
their bearing plates. Correction of these deficiencies is expected to lead 
to an accuracy much better than one percent. 

The error analysis of vertical c.g. determination (Appendix A) indi- 

2 

cates that the chosen pendulous method of determining is not a precise 
method of locating center of gravity. Nor should it be. By intent, the test 
parameters were chosen to place minimal dependence on knowledge or 
determination of c. g. position: 


dD 

dr 


0 


The minimum of 1T D, IT or of period of oscillation, is much more precisely 
bounded than is the value of n r M at which this minimum occurs. 
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Table 1. Helios prototype spacecraft mass property estimates 
provided by Messerschmitt- Bolkow - Blohm, GMBH 
(Tentative information) 


w = 

358 kg 


z = 

980 mm 

(above separation plane) 

I 

ox 

165 kg • 

2 

m 

I 

157 kg • 

2 

m 

oy 


I 

oz 

192 kg . 

2 

m 
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Table 2. Timing measurements for fixture plus adapter and V-band 


Hole 

Designation 

X 

mm 

Initial 

Amplitude 

mm 

Elapsed Time 
for 40 Cycles 
sec. 

Period, T, at 
'’Zero Amplitude” 
sec. 

B 

0 

25 

113. 739 




25 

113. 738 




13 

113. 727 




13 

113. 725 




0 

113. 720 

2. 8430 

B (repeat) 


25 

113. 724 




25 

113. 718 




13 

113. 706 




13 

113. 707 




0 

113. 700 

2. 8425 

c 

43. 18 

25 

113.608 




25 

113.605 




13 

113. 592 




13 

113. 593 




0 

113. 585 

2. 8396 

C (repeat) 


25 

113. 588 




25 

113. 582 




13 

113. 593 




13 

113. 592 




0 

113. 595 

2. 8399 

D 

86 . 36 

25 

113. 572 




25 

113. 570 




13 

113. 567 




13 

1 13. 562 




0 

113. 560 

2. 8390 

D (repeat) 


25 

113. 541 




25 

113. 540 




13 

113. 531 




13 

113. 534 




0 

113. 532 

2. 8383 

E 

128. 0 

25 

' 113.563 




25 

113. 566 




13 

113. 555 




13 

113. 556 




0 

113. 552 

2. 8388 
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Table 2 (contd) 


Hole x 

Designation mm 


Initi al 
Amplitude 
mm 


Elapsed Time 
for 40 Cycles 
sec. 


Period, T, at 
"Zero Amplitude" 
sec. 


E (repeat) 


F 166.9 


F (repeat) 


25 

1 13. 604 

25 

113. 600 

13 

113. 590 

13 

113. 591 

0 

113. 586 

25 

113.697 

25 

113.698 

13 

113.682 

13 

113.683 

0 

113.678 

25 

1 13. 722 

25 

113.718 

13 

113. 716 

13 

1 13. 714 

0 

113. 714 


2. 8397 


2. 8420 


2. 8429 


Initial estimates for statistical determinations: 


2 

^o 


1.035 m 2 


r = 0.925 m 


JPL Technical Memorandum 33-707 


23 



Table 3. Timing measurements for fixture plus spacecraft 
Oscillation about S/C x-x axes 


Hole 

Designation 

X 

mm 

Initi al 
Amplitude 
mm 

Elapsed Time 
sec. 

Period, T, at 
"Zero Amplitude' 1 
sec. 

D 

86. 36 

25 

106. 879 




25 

106. 878 




13 

106. 865 




13 

106.872 




0 

106.865 

2.6716 

D (repeat) 


25 

106. 860 




25 

106. 858 




13 

106.840 




13 

106. 843 




0 

106. 835 

2. 6709 

c 

43. 18 

25 

106.607 




25 

106. 594 




13 

106. 593 




13 

106. 595 




0 

106. 592 

2.6648 

C (repeat) 


25 

106. 544 




25 

106. 539 




13 

106. 528 




13 

106. 525 




0 

106. 522 

2. 6631 

B 

0 

25 

106. 376 




25 

106. 374 




13 

106. 340 




13 

106.342 




0 

106. 328 

2.6582 

B (repeat) 


25 

106.417 




25 

106. 425 




13 

106.403 




13 

106.400 




0 

106.387 

2. 6597 

B (repeat) 


25 

106.429 




25 

106.425 




13 

106.419 




13 

106.419 




0 

106.413 

2.6603 
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Table 3 (contd) 


Hole 

Designation 

__ Initial 

Amplitude 
mm r 

mm 

Elapsed Time 
sec. 

Period, T, at 
"Zero Amplitude" 
sec. 

A 

-43. 18 25 

106.440 



25 

106.434 



13 

106.408 



13 

106.414 



0 

106. 395 

2.6599 

A (repeat) 

25 

106. 392 



25 

106. 391 



13 

106.381 



13 

106.382 



0 

106.379 

2. 6595 

j 

-91.31 25 

106. 568 



25 

106. 566 



13 

106. 549 



13 

106. 547 



0 

106. 540 

2.6635 

J (repeat) 

25 

106. 509 



25 

106. 512 



13 

106. 489 



13 

106.489 



0 

106.476 

2.6619 

J (repeat) 

25 

106.531 



25 

106. 521 



13 

106.494 



13 

106.497 



0 

106.480 

2.6620 

Initial estimates for statistical determinations: 



p 2 = 0. 8010 m 2 

r o 

r = 0.905 
o 

m 
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Table 4. Timing measurements for fixture plus spacecraft 
oscillation about S/C y-y axes 


Hole 

De signation 

X 

mm 

Initial 

Amplitude 

mm 

Elapsed Time 
for 40 Cycles 
sec. 

Period, T, at 
"Zero Amplitude" 
sec. 

D 

86. 36 

25 

106.986 




25 

106.980 




13 

106.965 




13 

106.960 




0 

106. 955 

2.6739 

D (repeat) 


25 

106. 987 




25 

106. 982 




13 

106. 956 




13 

106.962 




0 

106. 944 

2.6736 

c 

43. 18 

25 

106. 762 




25 

106. 761 




13 

106. 743 




13 

106. 741 




0 

106. 727 

2.6682 

C (repeat) 


25 

106. 808 




25 

106. 809 




13 

106. 789 




13 

106. 785 




0 

106. 779 

2. 6695 

B 

0 

25 

106.636 




25 

106.644 




13 

106.627 




13 

106.625 




0 

106.620 

2. 66 55 

B (repeat) 


25 

106. 595 




25 

106. 586 




13 

106. 567 




13 

106. 570 




0 

106. 558 

2.6640 

B (repeat) 


25 

106.663 




25 

106.636 




13 

106.612 




13 

106.609 




0 

106. 589 

2.6647 
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Table 4 (contd) 


Hole 

Designation 

X 

mm 

Initial 

Amplitude 

mm 

Elapsed Time 
for 40 Cycles 
sec. 

Period, T, at 
"Zero Amplitude” 
sec. 

A 

-43. 18 

25 

106.678 




25 

106.676 




13 

106. 646 




13 

106. 645 




0 

106.635 

2.6659 

A (repeat) 


25 

106. 647 




25 

106.648 




13 

106.632 




13 

106.638 




0 

106.634 

2. 6659 

j 

-91. 31 

25 

106. 783 




25 

106.782 




13 

106. 765 




13 

106. 758 




0 

106. 753 

2. 6688 

J (repeat) 


25 

106. 849 




25 

106. 841 




13 

106.827 • 




13 

106. 829 




0 

106. 820 

2. 6705 

Initial estimates for 

statistical determinations: 



2 

0. 8079 m 2 

r = 0.9035 

m 


o 


o 
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Table 5. Timing measurements for fixture plus spacecraft 
oscillation about S/C p-p axes 


Hole 

Designation 

X 

mm 

Initial 

Amplitude 

mm 

Elapsed Time 
sec. 

Period, T, at 
"Zero Amplitude" 
sec. 

D 

86.36 

25 

107. 107 




25 

107. 102 




13 

107. 080 




13 

107. 082 




0 

107. 070 

2.6768 

D (rerun) 


25 

107. 089 




25 

107. 096 




13 

107. 089 




13 

107. 090 




0 

107. 090 

2. 6773 

c 

43. 18 

25 

106. 810 




25 

106. 804 




13 

106. 789 




13 

106,784 




0 

106.779 

2.6695 

C (repeat) 


25 

106.870 




25 

106. 861 




13 

106.848 




13 

106.848 




0 

106.840 

2. 6710 

B 

0 

25 

106.699 




25 

106.691 




13 

106.682 




13 

106.682 




0 

106.676 

2. 6669 

B (repeat) 


25 

106. 723 




25 

106. 724 




13 

106. 704 




13 

106.707 




0 

106.698 

2. 6675 

A 

-43. 18 

25 

106. 759 




25 

106. 760 




13 

106. 754 




13 

106. 751 




0 

106. 752 

2.6688 
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Table 5 (contd) 


Hole 

Designation 

Initial 

x Amplitude 

mm r 

mm 

Elapsed Time 
sec. 

Period, T, at 
"Zero Amplitude" 
sec. 


25 

106. 722 



25 

106. 712 



13 

106.708 



13 

106. 712 



0 

106. 707 

2.6677 

J 

-91.31 25 

106.864 



25 

106, 850 



13 

106. 831 



13 

106.830 



0 

106. 820 

2. 6705 

J (repeat) 

25 

106. 892 



25 

106. 884 



13 

106. 860 



13 

106. 864 



0 

106. 852 

2. 6713 

Initial estimates for statistical determinations: 



p 2 = 0. 8097 m 2 
r o 

r = 0. 9042 
o 

m 
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Table 6. Processed results of oscillation tests 



Fixture plus adapter 

and V-band 




Best Fit 

Based on 
Static- Test 

r o (m) 


0.9161 

0. 9042 

2 < 

Po (m ) 


1. 0296 

1. 0297 

A 2 / 2 

A p / p 
r O r o 


0. 17% 

0. 17% 

Ar Q (mm) best fit and static test 

11.9 


Fixture plus spacecraft 

Axis 


Best Fit 

Consistent C. G. 

X-X r (m) 

o ' 


0. 9131 

0. 9065 

2 / 2 \ 
P Q ) 


0. 7972 

0. 7976 

A 2 , 2 

A P 0 /P 0 


0. 16% 

0. 16% 

Y-Y r Q (m). 


0. 9016 

0. 9065 

2 / 2 ^ 
P Q (m ) 


0. 8040 

0. 8038 

2 ? 
Ap "/ n 
r o ' M o 


0. 19% 

0. 11% 

P-P r (m) 

o ' ’ 


0. 9049 

0. 9065 

2 / 2 . 
P Q ( m ) 


0. 8061 

0. 8060 

A 2 , 2 

A Po ; Po 


0. 19% 

0. 19% 

r Q (mm) 

average 

906. 5 


Ar^ (mm) 

maximum deviation 

6 . 6 


Ar^ (mm) 

standard deviation 

5. 2 
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Table 7. Timing measurements for bare fixture plus adapter 
and V-band oscillation about S/ C roll axis 
(initial amplitude, G q = 0. 016 rad. ) 


AI 2 
kg • m 

Elapsed Time 
for 80 Cycles 
sec. 

Average 
Period, T, 
sec. 

T 2 

2 

(sec. ) 

0 

144. 872 

144. 873 

1. 81090 

3. 2794 

1. 0192 

145. 818 
145. 813 

1. 82270 

3. 3223 

2. 0380 

146. 761 
146. 758 

1. 83450 

3. 3654 

3. 0575 

147. 684 
147.684 

1. 84605 

3.4079 

4. OJ65 

148.623 
148. 633 

1. 85784 

3. 4516 

5. 0961 

149. 555 
149. 562 

1. 86949 

3.4950 
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Table 8. Timing measurements for fixture plus spacecraft 
oscillation about S/C roll axis 
(initial amplitude, 6 = 0.016 rad.) 


AI 2 
kg • rn 

Elapsed Time 
for 80 Cycles 
sec. 

Average 
Period, T, 
sec. 

T 

(sec. ) 

0 

270.671 




270.666 

3. 38336 

11.4471 

4. 7407 

273. 013 




273.017 

3.42169 

11.6464 

9.4801 

275. 355 




275. 370 

3.44204 

11. 8476 

- 

275. 301 




275. 324 

3. 44141 

11. 8433 

14. 219 

277.644 




277.673 

3. 47070 

12. 0458 

18. 958 

279. 981 




279. 979 

3.49975 

12. 2483 


279. 938 




279. 933 

3.49920 

12. 2444 
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Table 9. Results of roll moment of inertia determinations 



Configur ation 

Iz 7 
(kg . m ) 

AI /I 
z z 

(%) 

1. 

Fixture + spacecraft 

271.49 

0.092 

2a. 

Fixture as measured 

77. 547 

0, 212 

2b. 

Increment for fixture with spacecraft# 

0. 045 


3. 

Spacecraft 

193.90 . 

0. 154 

^Longer bolts were used for attaching weight increments. 
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Fig. 3. Pendulous support for 
spacecraft 


Fig. 4. Cradle ballast weight with 
optional positions for pivot pin 
(symmetrical about cradle 
centerplane) 
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Fig. 5. Spacecraft on cradle with 
pins in hole ,? J M (see Fig. 8) 


Fig. 6. Frontal view of spacecraft 
on cradle 
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Fig. 7. Photoelectric sensor in 
place under cradle 


Fig. 8. Electronic 
counter /timer units 
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Fig. 9. Pivot-pin hole geometry relative 
to spacecraft separation plane 
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r, mm 


Fig. 11. Test scatter about statistically 
derived values of 

r o 
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Fig. 12. Roll moment-of -inertia 
fixture with adapter and V -band 
installed 


Fig. 13. Lower flexure -pivot 
assembly as viewed from above 
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Fig. 14. Installation view of 
spacecraft on roll moment-of- 
inertia fixture 
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Fig. 16. Installation of incremental test 
weights on transverse arm 
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Fig. 18, Plot of T2 vs AI for fixture 
plus spacecraft 
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APPENDIX A 


ERROR ANALYSES RELATED TO DEFINITION OF THE 
INERTIA ELLIPSE IN THE x-y PLANE 


Weights 

In the weighing operations involving fixturing, either of two load cells 
was used. These load cells were of 454-kg and 908-kg ratings, respec- 
tively. They were calibrated on a certified deadweight tester located in 
one of JPL's instrument laboratories to a read-out accuracy of 0. 1 lb. , 
or 0. 046 kg . Let 


Then 


W— = weight of the cradle, plus adapter and V-band 
r 

W, - weight of the hoisting sling 
h 

W = combined weight, 
c 


W, 


w - w, 

c h 


(A. 1) 


dWj 

WZ 


i 


'9W t 


w „ l aw 


dW + w 
c W 




aWj 

aw7 


dW, 


(A. 2) 


or 


AW, 


W, 


/W c \AW c /W h \ 
l W F.i W c \ W f/ 


AW } 


max 


From the calibration data, the following estimates are made: 


(A. 3) 


AW 

c 

W 


= 0 . 0002 ; 


AW 

~w7 


h 


= 0.001 
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Thus 


AW, 


W, 


max 


0. 00035 or 0. 035% 


Now, let 


W T = total weight of fixture plus spacecraft plus 
hoisting sling 

Wg = weight of spacecraft 


w s = w T - w F - W h 


(A. 4} 


By treatment similar to that accorded W_, it can be shown that 

r 



= 0.0010 or 0.1% 


(A. 5) 
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Static Determination of Fixture C.G . 

In addition to errors in forces and distances, errors are also assumed 
in angular alignments, G and 4>> as depicted in the sketch below. 


I / p 



P = W r cos 6 - aG 

cos 4> L cos G + b sinG ’ 

P 3 A 

r = - cos 4> (L + b tan G ) - (A. 7) 


dr _ P /_9r_\ dP W (dr\ dW 
r r iap/ P + r \aw/ W 


L / 9r \ dL b / 9r \ db _a_ / 9r \ da 
r \3L/ L. r \9b/ b + r \8a/ a 
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9r 

9P 


COS 4> 

= W 


( L + b tan G ) 


9r _ P cos 4 > , T , , . Q . 
aw w 2 (L + b tan 9 ) 


3r P , 

"aTT = -w cos<|) 


3r 

ab 


cos 4 * tan 0 


dr _ e 

9a ~ cos G 


0 r 

ae 


Pb , 2_ a 

- 777 - cos d> sec G ~ 

W cos 0 


(1 - G tan 6 ) 


8 r 

a4> 


p 

-yjr sin 4 ) |L + b tan G 


Here, 

P = 186. 71 kg (less V-band) 

W = 417. 72 - 7. 54 = 410. 18 kg (less V-band) 
L = 2. 169 m 
b = 0. 583 m 
G = ±0° 3 1 (0. 00087 rad. ) 

4 > = ±0° 30’ (0. 0087 rad. ) 
a = 0. 0159 m 
r = 0. 9041 m 
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Then 


Ar 

r 


1. 09 



1. 09 


AW 

W 


+ 1.09 


AL 

L 


with 


AP 

P 


0. 0007; 


AW 

W 


0. 00035; 


AL 

L 


0 . 0006 ; 


— = 0.0018, or 0. 18%; Ar = ±4 mm 
r 


"Least Squares" Determination of the Radii of Gyration and Center of 
Gravity for the x-x, y-y and p-p Axes 

From Eq. (1.15), the period T is given by 


T 



where a and g are known and the distance, r, from pin support to the 
system c. g. is given by 


r = r + x 
o 

r is a reference length here, equal to the distance from the base of 
"Pin B" to the center of gravity, and x corresponds to increments on r Q 
to the different pin holes on the fixture. 

The distance, r, from pin support to the c. g. was varied and the cor- 

2 

responding periods noted. The radius of gyration p Q is first determined 
from the minimum value of T in a plot of T versus x. The following pro- 
cedure determines the center of gravity position r and radius of gyration 
2 

p for the best fit of the experimental data. 
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To this end, the mean square error 


2 

<r 





2 


where 


r . = r + x. 

1 o i 


is minimized with respect to and . 

In order to simplify the algebra, the following was considered: 



r . = r + x. 

1 o 1 


with 


4TT 2 /g 


D 


Then 


D. 

i 



r . + a 
i 


D. + 2a 
i 


(r + a + x. ) + 
o i r 


2 ^ 

p + a 
r o 

+ a + x. 
o 1 
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Put 


R = r + a; 
o 



E. = D. + 2a 
1 1 


Then 


E. = 

i 


R + X i + R 


B 


+ x. 


and the mean square error 


/ 


\ 

n ( 



2 1 ) 

, B 

1 

' ■ Ei ' 

{ R + x, + -p 
v x R + x. 


( 

i 



is now considered. 

The minimization with respect to R and B results in two nonlinear, 
simultaneous equations which are very difficult to solve. The problem is 
linearized by considering small increments on initial values (very close 
approximations) for R and B. Bet 

R = R + a; B = B + (3 

o o 


and 


R = R q (1 +i), B = R q 2 (B q + n ) 


where R and B 
Q O 


are the close approximations to R and B and 


r, = «/R o , 



B 

o 


B /R 2 
o o 
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then 


R -f a + x, 
o 1 


B + p 

Q 

R + a + x. 
o i 


(A. 9) 


° M L J 1 + l + * t 


where 


E. = E./R ; x. = x. /R 

i 1 o 1 i o 


<r 1 )- - ( B 0 + n ) / r 

R- - nE E i - 


(1 + X, 


1 + X. 
1 


Then, neglecting terms of order t/ 


ir-*Ek 


(1 + x.) - 


(i + x. ; 


(1 + X.) 


(i + x,: 


(A. 10) 
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or 


where 


n 


<r j.V' 
R “ nZj 


(a. + b.f, + c.ri) 
1 l 1 1 


i= 1 


(A. 11} 


a = E - (1 + x. ) ^ = - 

1 (1+x.) 


B 

o 

(1 +X.) 


1 


c. 

1 


-1 

(1+X.) 


Differentiating Eq. (A. 11) with respect to f, 


and since 




2b (a. + b. 4 + c.ri ) 
i i m l 1 


d 

a; 



o 


for a minimum 



n 

n 

n 


^L b i 2 + 


c. + / b. a. = 0 

1 jLmJ 1 1 


i=l 

i=l 

i= 1 


(A. 12) 
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Similarly for r\ : 


and 


n 


a 

^ U, 


= — / ^ 2c. (a. b.t, + c.n ) 
n / 4 11 l 1 


i= 1 


9t 1 


o~ 

R. 


= 0 


n 


n n 




i= 1 


i = 1 i=l 


(A. 13) 


Solving for t, and t| ; from (A. 12), 


n = 


^h" (Eh *, + i XV 


{A. 14) 


and from Eqs. (A. 13) and (A. 14) 


£ = 


y* a. b, \ c. 2 - y' a. c. y'b. c. 
L-j i i Z-* i i i Z-/ i i 

(Xh =i) 2 - X'V 2 a 


(A. 15) 


Thus, knowing the values of t, and r| , R and B are found from the 
relationships 


R = R + a and B = B +3 

o o 


This was further developed into the following iterative scheme suitable 
for a small computer program. 
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Values of £ and r\ are calculated for the close approximations to R 
and B, and respectively, and the values of R and B calculated there- 
from are substituted into the expressions for £ and p as new close approxi- 
mations. The new values of ^ and p are then used to find new values of R 
and B and the cycle repeated until sufficient convergence is achieved. 

That is: 


R ^ = R (j) + a (j) 


B (j + - R (j) + 0^) 

n 


A small computer program was written. Convergence to numerically 

2 

acceptable values of r^ and was realized typically within five to six 

iterations . 


Radii of Gyration with Consistent C. G, Location 

An additional step was taken following the determination of the center 

of gravity location and radius of gyration of each of the three axes tested 

(i. e, , x-x axis, y-y axis and p-p axis of the spacecraft and fixture). This 

step makes use of the fact that the center of gravity position should be the 

same for all three axes. The values for r were close but not exactly the 

o 

same for the three axes and, accordingly, an average value of the three 

values of r was taken to give a consistent value of r , This value of r was 

° 2 ° ° 
then used to compute new values of p Q for each of the three axes as follows: 

from before, 

2 2 

? ? p + r 

D = (g/4ir )T and D = — ; ; r = r +x 

' 6 ' r + a o 


Define the radius of gyration p Q for each location of the axis of rotation: 


p . ^ = (r. + a) D. - r . ^ 

l i i l 


(A. 16) 


i = 1 , • • • , n 
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Then 


p is the mean of the n values of p. 
r o r i 


n 


n 




i=l 


i = l 


(A. 17) 


which can be written as Eq. (1. 21). 

2 2 

Finally, the variance S on p is calculated by 


— ~r fp 2 - (r. + a) D. + r.‘ 

n - 1 l r o v l ' l i 

i = l X 


(A. 18) 


for the best fit. 

A similar technique was used in the case of the bare fixture where the 
center of gravity value calculated from the static test was used in preference 
to the one determined by experiment and best fit. In all cases the differences 
in the radii of gyration calculated by the best fit program and the consistent 
C. G. location technique were small (0. 03% in the worst case). The consis- 
tent C. G. location technique using Eq. (1.21) is believed to give the best 
determination of the radii of gyration. 

Extraction of Spacecraft Inertial Properties 

The variances on the radii of gyration of the spacecraft alone were then 
calculated from: 


S s 2 = (1 + M p / Mg) S T 2 + M F /M g S F 2 
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where 


= mass of fixture 
r 

Mg = mass of spacecraft 

= variance on fixture radius of gyration 

= variance on fixture + spacecraft radius of gyration 

= variance on spacecraft radius of gyration 

The error estimates on the spacecraft moments of inertia are then 
given by 



I 


S 


M, 


2 

S 


and 


or 



AM S Ap o 2 



o 



(A. 20) 


N. B. Alg/Ig is an estimate of the error. Alg/Ig is made up of an absolute 
error on Mg and a statistical error estimate on p^g and therefore not 
strictly the maximum absolute error. 
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Axis 


P 0 2 (M 2 ) _S S (M 2 ) 


x-x 

0 . 

5201 

0 . 

0026 

y-y 

0 . 

5341 

0 . 

0022 

p-p 

0 . 

5389 

0 . 

0031 

Finally, 

from Eq. 

(1. 25) 




I (kg • m 2 ) 

181. 7 
186,6 
188. 3 


r 


oS 



r oS' + r oT 


(%) 

0. 61 
0. 52 
0. 67 


=> r QS = 0.9093 m and Z = 1838.45 mm 

The Inertia Ellipse 

From Eq. (1.29) 


where 



au 

"ai 




y 


tan a 


0. 7483 


au 

~di 




P 


1 

sin 2 a 


-1. 082 
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au 

-rr-^L = (i - I ) - 2U ctn 2 a = 0 
3 a ' y x xy 


Thus 


AU 


U 


J£L - 


= 0.3340 


xy 




+ 0. 7483 


te) 1 ' 


AI 

y 


1. 082 


/ 1 \ AI 

ITT 1 

\ xy/ p 


= -0. 1072 - 0. 2102 + 0. 3951 


= 0.0777 


For the error in <£ , 


d<t> 

4> 


= ♦ («%-) dU *y + * ) 


dl 


x 


1 / 54 > \ 

+ KJ 


dl 


F rom Eq . (1.31) 


9<t> 


- (I - I ) 
x y 


3U , T T . ATT 2 

xy (I - I ) + 4U 

x y xy 


34 * _ _ 94 * _ 


u 


31 


31 


J£L 


y (i 


2 2 
I ) + 4U 

y xy 


(A. 22) 
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Thus 


= 0.052 
9 


or 


A4> = ±1.4' 


Finally, 


dI r 1 / 9I r\ 1 / 9I r\ 

i i (ai ) dI x + i ( 9 1 ) 

r r \ x/ r \ y/ 


dl 


WJM 

^ \ 8U J 


dU xy + -T (-5?-) d * 


(A. 23) 


From Eq. { 1 . 28) 


91 

r 

- cos <(> ; 

CP 
i— i 

>■* 

91 

X 

di 

y 

91 

r 

= sin 2<j>; 

91 

r 

9U 

xy 

9«|> 

AI 


/' 

AI 


• 2 a 

- sin <p 


rs 


T ^1 

-j = 0.6 5%; — J - = 0.47% 


Errors in Spacecraft C. G. Determination 
From Eq. (1. 25), 


M, 


tt) + r . 


oS M g ' oT oF oT 


dr 


9 r q 9r ~ 9r 9r 

oS ” 3M„ dM F + 3M, dM S + 9r ^ dr oT + 8r °„ dr oF 


oT 


oF 
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Ar 


oS 


M f AM f M f AMg 

= M (r oT ‘ r oF^ M + Ml (r oT " I ’oF ) 
max o Fob 


/M \ M 

+ (^7 + ') Ar oT + Ar oF 


Here, 



r 


oT 


r 


oF 



AM S 

M s 


Ar 

Ar 


oT 

oF 


and 


4 16 . 80 kg 
350. 08 kg 

906.5 mm (average value from Table 6) 
904. 2 mm (from static test) 

0. 0035 
0 . 001 

6,6 mm (max. deviation from Table 6) 

4.4 mm (from error analysis of static test) 


Ar _ = ±19. 3 mm 
oS 


(A. 24) 
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APPENDIX B 


ROLL MOMENT OF INERTIA DETERMINATION 
AND ERROR ANALYSIS 

Roll Moment of Inertia Determination 

From Eqs. (2.3 and (2.4) we have 

T 2 = T q 2 + mAI (B. 1) 

where 

T 2 

m = ~Y~ (B. 2) 

o 

The period, T , for zero mass increment and the slope m can be 
determined statistically from a least squares fit of the test data to a plot 
of T versus AI. 

The mean square error is given by 


2 

cr 



Minimizing with respect to m and 


T 

o 


2 


3cr 


3m 


= 0 




mAI. 



(B. 3) 


Then, 


T 


2 


o 



AI. 

l 


+ m 



(AI.) 2 



T. 

l 


2 


(B. 4) 
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and 



0 



mAI. (-1) 


which gives 

t 0 2 n + m E Ai i = S T i 2 


Dividing Eqs. (B*4) and (B. 5) by m: 



From Eq. (B. 6) , 



and substituting Eq. (B„ 8) into Eq. (B.7) it is found that 


m 


= I 


S T i' SK ) 2 - E AI i EV 
n £ AL Tj 2 - X T i 2 E 


AI. 


Error Analysis 


2 . 2 

The variance <r on the period squared T is given by: 


2 

cr 



T 


2 

measur 



2 


(B. 5) 


( B. 6 ) 


(B.7) 


(B. 8) 


(B. 9) 


(B. 10) 
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2 2 

The slope and intercept variances, o - and <x £ respectively, 
then follow: ^o 


m 


2 


EK - AI ) 


(B. 11) 


2 

a 2 
T 

o 


E K ) 2 

n £ K - AI) 2 


(B. 12) 


The variance on the roll moment of inertia I Q is given by: 


m 


ini = fa T - fa m 
o o 


AI 

o 


AT 2 

o 

1 

Am 

I 

o 


T 2 
o 

T 

m 


or 


AI 


AT 


o Am 
m 


max 


t °“ T 2 

I T cr 

o o + m 


X o T 2 m 
o 


(B. 13) 


The error estimate on the roll moment of inertia of the spacecraft 
alone is then given by 


2 2,2 

<r = <r + <r 

oS oT oF 


(B. 14) 
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where subscripts are 
S - spacecraft 
T = spacecraft + fixture 
F = fixture alone. 
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APPENDIX C 


MODEL TESTS AND ASSOCIATED ANALYSES 

The objectives of the model tests were twofold: 

(1) to acquire knowledge that could reflect to advantage in the 
design of the full-scale test fixturing. 

(2) to assess the accuracy potential prior to start of fabrication 
of the full-scale fixturing. 

The implementation of this model program was predicated on the 
assumption that the moments of inertia of prismatic bars representing 
"fixture" and "spacecraft" could be calculated with sufficient accuracy to 
serve as "absolute standards" for judging experimental accuracy. 

Accordingly, bars of precision-ground tool steel were machined and 

drilled according to design requirements, and actual dimensions were 

checked to an accuracy of 0.003 mm by JPL's Inspection Department. The 

o 

actual dimensions were used in calculations of p of each element with 

o 

allowance for material removed for containing the pivot pin. 

Figure C-l is a view of the "fixture" simulator in place on its sup- 
porting structure, which provided appropriately spaced drill rods as the 
ways on which the pivot pin rode. 

Figure C-2 shows the elements of the complete model; the combination 
of the two smaller bars represents the "spacecraft." 

The test models were purposely designed so that they could be used 
as "reversible pendulums," and oscillation tests were conducted in both 
aspects. Actual weights of the model elements were measured on a certi- 
fied gram balance. 

The test data were processed in the same manner as that accorded the 
data of the full-scale tests. The results are summarized in Table C-l. 

The experimentally determined value of the "fixture" I is in excellent 

o 

agreement with the "actual" value, which was derived from a measured 

2 

weight and a value of p Q calculated from the bar geometry. The corre- 
sponding comparison for "fixture plus spacecraft" is not nearly as good, 
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suggesting some systematic inaccuracy that was not revealed by the statis- 
tical error analysis. 

It is to be noted that, under the combined weight of "fixture" and 

"spacecraft, " galling of the 4. 8-mm (0. 188-inch) diameter drill rod that 

served as the pivot pin was observed. This situation appeared to be relieved 

by using a plug gage of greater surface hardness. However, it is not certain 

that the problem was completely eliminated, as the bearing stresses under 

nominal "point contact" were very high. Any flattening of the pivot pin 

increases the pin radius, a, reduces the parameter, D , and hence the 

Z m 

experimentally derived value of p . Thus the design of the full-scale 
fixturing provided for line contact of the pins on flat ways, although this 
ideal situation was not completely realized in the full-scale tests for reasons 
noted. 

The error in the "spacecraft" moment of inertia shown in Table C-l 
is, of course, magnified in the process of removing the "fixture" from "fix- 
ture plus spacecraft, " Even so, an accuracy of better than 1% was achieved, 
demonstrating the high accuracy attainable by the method. 
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Table C-l. Comparisons of measured and ’’actual 11 
inertial properties 


1. ’’Fixture" 

Mass (kg) 

2 

Experimentally determined I (kg m ) 
Estimated statistical test error 

f 2 

"Actual” moment of inertia, 1^ (kg m ) 

r 

Experimental error (%) based on I 

o 

2. "Fixture 11 plus "spacecraft" 

Mass (kg) 

2 

Experimentally determined I q (kg m ) 
Estimated statistical test error 
"Actual" moment of inertia, I q (kg m ) 

i 

Experimental error (%) based on I 

3. "Spacecraft" 

Mass (kg) 

2 

Experimentally determined I (kg m ) 
Estimated statistical test error 

i 2 

"Actual" moment of inertia, I (kg m ) 

I 

Experimental error (%) based on I 


4. 622 
0. 08456 
±0. 23% 

0. 08458 
-0. 03% 


8. 587 
0. 13270 
±0. 24% 

0. 13324 
-0. 40% 


3. 965 
0. 04814 
±0.43% 

0. 04852 
-0. 78% - 
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Fig. C-l. Model for test of radius- 
of-gyration measurement concept 



Fig. C-2. Elements of 
complete model 
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